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Abstract — Superconducting metamaterials combine the advan- 
tages of low-loss, large inductance (with the addition of kinetic 
inductance), and extreme tunability compared to their normal 
metal counterparts. Therefore, they allow realization of compact 
designs operating at low frequencies. We have recently developed 
radio frequency (RF) metamaterials with a high loaded quality 
factor and an electrical size as small as ~ A/658, (A is the free 
space wavelength) by using Nb thin films. The RF metamaterial 
is composed of truly planar spirals patterned with lithographic 
techniques. Linear transmission characteristics of these meta- 
materials show robust Lorentzian resonant peaks in the sub- 
100 MHz frequency range below the T c of Nb. Though Nb is 
a non-magnetic material, the circulating currents in the spirals 
generated by RF signals produce a strong magnetic response, 
which can be tuned sensitively either by temperature or magnetic 
field thanks to the superconducting nature of the design. We have 
also observed strong nonlinearity and meta-stable jumps in the 
transmission data with increasing RF input power until the Nb 
is driven into the normal state. We discuss the factors modifying 
the induced magnetic response from single and 1-D arrays of 
spirals in the light of numerical simulations. 

Index Terms — RF metamaterials, superconducting devices, ar- 
tificial magnetic response. 



I. Introduction 

METAMATERIALS are engineered structures composed 
of sub-wavelength artificial unit cells that form a homo- 
geneous effective medium for the electromagnetic waves prop- 
agating through them. This effective response can be described 
by a coarse-grained electric permittivity and the magnetic 
permeability Q. By manipulating these two electromagnetic 
parameters, interesting properties unavailable in nature have 
been demonstrated, such as negative refractive index (2) and 
artificial magnetism 131, |4|. 

Most metamaterials are designed to give a magnetic re- 
sponse at microwave or higher frequency bands O, Q. 
They usually employ split ring resonators (SRRs), or their 
derivatives, and can be generally understood with a lumped- 
element inductor-capacitor circuit model 0. 
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As the wavelength, A, of the electromagnetic waves in- 
creases towards the radio frequency (RF) band, miniaturizing 
the inclusions forming a metamaterial becomes more challeng- 
ing (6|. The Swiss Roll metamaterial |6], a three dimensional 
(3-D) geometry which inspired the SRR design, has achieved a 
reasonable magnetic response in the sub- 100 MHz band. Since 
the design is implemented with normal metals (thick Cu sheets 
separated by dielectrics), ohmic losses and bulkiness are the 
limiting factors. 

The most important usage of RF metamaterials reported so 
far is in Magnetic Resonance Imaging devices Q; however 
they can function in a wide variety of applications including 
magneto-inductive lenses (8), RF antennas |9], filters and 
compact resonators |[T0li . 

II. The Design 

Our Nb RF metamaterials are constructed from micro- 
structured 2-D spirals. Earlier experiments ifTTIl . fT2ll utilized 
planar spirals made of thick Cu films (~35 |Ltm thick in IfTTIl 
and ~0.25 mm thick in fT2ll ) on dielectric substrates. Such 
thick coatings are required to minimize Ohmic losses, making 
the design intermediate between 2-D and 3-D. Further reduc- 
tion of the resonant frequency of such spirals is not possible 
due to the scaling of Ohmic dissipation with spiral length. 

Because the electrical size of the elemental spiral is signif- 
icantly smaller than the operating wavelength, a quasi-static 
approach can be used, and the design can be based on a 
lumped element analogy as in SRRs 1 1 3 ] . The large number 
of turns gives a significant geometrical inductance, L g iTUfl] 
and the spacing between them provides substantial capaci- 
tance, C 021 • Therefore, the geometry can be kept compact 
enough to achieve a sub- 100 MHz resonance frequency, fo = 
l/(27Ty/ (L g + Lk)C (despite the large A). Since the design 
is realized with superconducting films, there is an additional 
term called kinetic inductance, 1 16] which is closely linked 
to the superfluid density. For example, when the transition 
temperature, T c , is approached, the superfluid density goes to 
zero, and can make the dominant contribution to the total 
inductance, resulting in a significant reduction in fo. 

III. Sample Preparation 

The fabrication starts with the deposition of 200 nm Nb thin 
films on 350 |Ltm thick quartz substrates, 3" in diameter, by 
the RF sputtering technique |17|. The transition temperature 
of the film is found to be T c ~ 9.25 K from resistance 
measurements taken by a Physical Property Measurement 
System (PPMS) |[T8ll . The main panel of Fig. 1 shows the 
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Fig. 1. (Color online) Main panel: Resistance vs. temperature of the 
as-deposited Nb thin film. Top-left inset: Zoom in on the resistance vs. 
temperature curve showing the transition in detail. Bottom-right inset: Optical 
photograph of a representative Nb spiral with D Q = 6 mm showing 40 turns 
in detail. 



resistance of the Nb film for a wide range of temperature 
from 300 K down to 1.5 K, whereas the top-left inset focuses 
on the transition region. The Nb film is patterned into a spiral 
geometry by photolithography and reactive ion etching. The 
etching is performed by a chemically active mixture of CF 4 /0 2 
(10% 2 ). The wafer is diced into chips having single spirals 
and into 1-D arrays composed of various numbers of spirals. 
An optical image of one of the single spirals giving sub- 100 
MHz resonant features with an outer diameter D Q of 6 mm 
and N = 40 turns is shown in the bottom-right inset of Fig. 1. 
The width w of the stripes and spacing s between them are 
equal and each is 10 |Ltm. 

IV. Building Blocks of the RF Metamaterial: The 
Single Spiral 

A. Transmission Data 

The RF transmission measurements are carried out on single 
elements of our RF metamaterials in an evacuated cryostat 
at temperatures ranging from 300 K down to 4.2 K. The 
sample is mounted between two RF magnetic loop probes, 6 
mm in diameter, whose axes are parallel to each other and 
the Nb spiral, as shown in the inset of Fig. 3. The spiral 
resonator is excited by the top probe providing the RF signal 
and the transmission is picked up by the bottom probe. Both 
loop probes are connected to an Agilent E5062A RF network 
analyzer. 

The full transmission spectrum on a single spiral with D Q 
= 6 mm, N = 40, w = s = 10 \xm is shown in Fig. 2(a). 
The data are taken at T = 4.3 K with an input power of 10 
dBm and demonstrate a fundamental resonance at 74 MHz as 
well as harmonics located above 100 MHz. The RF current 
density in the spiral at those resonance modes are examined 
with laser scanning microscopy (LSM) 1191 , |20l in Fig. 2(b)- 
(f). These images can be interpreted as showing approximately 
the local current density squared flowing in the windings of 
the spiral (light is large, dark is small current). The inner and 
outer borders of the spiral are marked by dashed circles. In 
the fundamental mode, an intense current distribution flows 
through the middle windings, and gets weaker near the edges 
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Fig. 2. (Color online) (a) Measured transmission |S2i| vs. frequency on a 
single spiral of D Q = 6 mm, TV = 40, w = s = 10 \im at T = 4.3 K showing the 
fundamental and higher mode resonances, (b)-(f) Laser Scanning Microscopy 
images on the same spiral showing RF current distributions inside the turns 
corresponding to the resonant modes in the transmission data. 



of the spiral (Fig. 2(b)). For the 2 nd harmonic (Fig. 2(c)), there 
are two strong current distributions, presumably of opposite 
sign, resulting in a partial cancellation of the fields, and 
responsible for the smaller transmission amplitude seen in 
Fig. 2(a). Higher harmonics (Fig. 2(d)-(f)) show an increasing 
number of large- amplitude circles, suggesting that the spiral 
acts almost like a distributed circuit resonator carrying integer 
multiple half-wavelengths of current in each standing wave 
eigenmode. 

The transmission data of the same spiral at the fundamental 
mode for two different loop- sample distances, h, are shown in 
the main panel of Fig. 3. Both data sets are taken well below 
the T c of Nb; data with h = 3.25 mm (filled stars) obtained at 
4.8 K and with h = 13.5 mm (empty circles) at 4.4 K. Since 
Nb is in the superconducting state, and the driving power is 
low (P = -10 dBm for the data shown in Fig. 3), the resonance 
peaks are independent of RF input power and are of Lorentzian 
shape (the solid curves are fits to the data). 

As seen from the figure, larger separation between magnetic 
loops leads to sharper transmission peaks, and consequently to 
a higher loaded quality factor, Ql, which is calculated as ifTTTl 
Ql = fo/5f(i.5dB)> where fo is the frequency at which the 
peak is located (i.e. resonance frequency) and 5f(\.$dB) is the 
bandwidth where IS21I is 1.5 dB below its peak value (see 
the top-left inset of Fig. 5). The loaded quality factor of the 
transmission data taken for h = 3.25 mm is ~1502 and for 
h = 13.5 mm one finds Ql ~7800. These quality factor values 
are up to several orders of magnitude larger than their normal 
metal counterparts fTTl . This brings up a discussion of the 
effect of coupling between the RF probes on the transmission 
| £21 1 data. 

B. Coupling Between RF magnetic Loops 

The measured Ql is determined by external and internal 
loss factors and can be written as Ql = Uuj/P t where U is 
the stored energy in the spiral on resonance, uj is the resonant 
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Fig. 3. (Color online) Main panel: Measured transmission IS21 1 vs. frequency 
on a single spiral of D = 6 mm for two different probe distances to the 
sample, h. Filled stars are data obtained for h = 3.25 mm, and empty circles 
are those for h = 13.5 mm. Solid curves are Lorentzian fits, which agree 
closely with the data, indicating linear response in both experimental data sets. 
Inset: Sketch of the arrangement for the microwave transmission experiments; 
a single Nb spiral sandwiched between two magnetic loops extending from 
coaxial cables. 

angular frequency and P t is the total dissipated power |23l . P t 
is the sum of power dissipated in the cavity (Ohmic losses), 
radiated power, and power leaking out of the spiral through 
the coupling loops. Since the measurements are done well 
below T c , ohmic and dielectric losses should be the smallest 
sources of dissipation. We find that the largest dissipation 
is due to the coupling losses from the RF loop probes. 
To confirm this, we have simulated the transmission curves 
numerically on a lossless single spiral of D Q = 6 mm having 
the same dimensions as those used in the experiments shown 
in Fig. 3, utilizing the High Frequency Structure Simulator 
(HFSS) l22l . As seen in Fig. 4, fo shifts to lower frequencies 
with increasing h, which is consistent with the data in Fig. 3, 
and the expectation from electromagnetic perturbation theory 
for a dominantly magnetic perturbation [24]. 

The loaded quality factor is observed to increase substan- 
tially as h increases. This is consistent with an increase in 
the coupling quality factor, Q coupling with weaker coupling, 
and Ql 1 = Q' 1 + Q coupling- A sigmoidal trend of Q L is 
shown in Fig. 5 along with the polynomial decrease in / 
with increasing h. 

V. Towards Creation of an Artificial RF Medium: 
1-D Arrays of the Spirals 

A. Linear Transmission Data below and above T c 

So far, the discussion has been focused on a single element 
of the artificial medium, namely a single spiral showing a 
sub- 100 MHz resonance. The metamaterial properties should 
be confirmed by examining an ordered collection of these 
elements. 

The excitation of the linear array is as follows; the top loop 
is aligned with the first spiral of the array whereas the bottom 
loop is aligned with the last spiral. The distance between the 
edges of two spirals in the array is 1.5 mm, and h is 3.25 
mm. The measured transmission spectrum on a 1-D array of 
7 spirals with D Q = 6 mm, N = 40, w = s = 10 \xm is shown 
in the main panel of Fig. 6. The Nb spiral array is excited 
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Fig. 4. (Color online) Calculated transmission IS21I vs. frequency by HFSS 
for a set of distances h between a single spiral of D Q =6 mm and one of 
the RF magnetic loops. As the driving loop is moved away from the spiral, 
fo shifts to smaller values and Q l becomes larger. 
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Fig. 5. (Color online) Main panel: Quality factor vs. driving loop probe 
distance to the spiral from HFSS calculations. The data (filled circles) follow 
a Boltzmann fit (solid curve). Top-left inset: The calculation of Ql for data 
at h = 6 mm (filled diamonds). Bottom-right inset: Resonance frequency 
decreasing with increasing probe distance to the spiral (filled stars). 

by an input power of dBm at T = 4.5 K, consequently 
multiple resonant peaks are observed which dissapear above 
the superconducting phase transition at T c . The inset of Fig. 
6 shows the calculated transmission spectrum on the same 
spiral array, and gives resonant features at a slightly higher 
frequency (HFSS does not include Lj~ or internal losses in the 
calculation). 

B. RF Power Evolution of Transmission Data and Nonlinear- 
ity 

In the linear regime, the behavior of transmission data 
does not show any significant changes with respect to RF 
power put into the spiral cavity. The superconducting phase is 
strong and RF oscillations are maintained without significant 
Ohmic losses (the spiral can be modeled as an ideal inductor- 
capacitor (L-C) resonant circuit). As RF power is ramped up, 
the magnetic field starts to penetrate the spiral, particularly at 
weak spots in the superconducting film [19]. For higher RF 
input power, the dissipation caused by both vortex dynamics 
and Joule heating degrades superconductivity and introduces 
a resistance R into the lumped element circuit analogue, 
and leads to nonlinear resonant features in the transmission 
data (251. 
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Fig. 6. (Color online) Main panel: Measured transmission IS21 1 vs. frequency 
spectrum on an array of 7 spirals with D Q = 6 mm, N = 40, w = s = 10 urn 
below and above the T c of Nb with an input RF power of P = dBm. Inset: 
The calculated tranmission on the same array by HFSS. 
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Fig. 7. (Color online) Measured transmission IS21I vs. frequency on the 
same array shown in Fig. 6 at various RF input powers ranging from P = -1 
dBm to P = 24 dBm. Increasing RF input power leads to non-linear features 
and unstable jumps in the transmission. Because an amplifier was used for 
these measurements, the background level increased compared to that seen in 
Fig. 6. 

Fig. 7 shows the measured transmission data on the same 
1-D array in Fig. 6 for a set of RF input powers, P at a 
fixed ambient temperature of ~4.2 K. At P = -1 dBm, the 
resonance peaks are sharp and have high Ql (linear case). 
With increasing input power dissipative processes nucleate at 
discrete sites in the spiral and increase the resistance R of 
the circuit, which yields nonlinear resonant features in the 
transmission spectrum. 

VI. Conclusion 

The linear and non-linear transmission characteristics of 
compact RF metamaterials employing Nb thin films are re- 
ported. In these RF metamaterials, the ohmic and dielectric 
losses are minimized by carrying out the experiments below 
the T c of Nb. The 2-D design simplifies the fabrication and 
maintains the uniformity of the artificial constituent elements. 
The robust magnetic response of the spirals can be sensitively 
tuned by manipulating the superconducting phase either with 
temperature or RF magnetic field. The resonant features in the 
spectra have an extremely high Q l (compared to those of al- 



ternative normal metal metamaterials), making them promising 
for RF applications. 
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